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The particle emission at intermediate velocities in mass 
asymmetric reactions is studied within the framework of clas- 
sical molecular dynamics. Two reactions in the Fermi en- 
ergy domain were modeled , ^*Ni-fC and ^*Ni-|-Au at 34.5 
MeV/nucleon. The availability of microscopic correlations at 
all times allowed a detailed study of the fragment formation 
process. Special attention was paid to the physical origin of 
fragments and emission timescales, which allowed us to disen- 
tangle the different processes involved in the mid-rapidity par- 
ticle production. Consequently, a clear distinction between a 
prompt pre-equilibrium emission and a delayed aligned asym- 
metric breakup of the heavier partner of the reaction was 
achieved. 

PACS number(s): 25.70 -z, 25.70.Mn, 25.70.Pq, 
02.70.Ns 



I. INTRODUCTION 

Heavy ion reactions in the Fermi energy domain (20 - 
100 MeV/nucleon) has been an intense field of research 
during the last decade. In particular, a lot of attention 
was paid to the strong emission of intermediate mass frag- 
ments (IMF) observed in the mid-rapidity range. [p|-p^. 
Lying in rapidity space between the contributions of the 
projectile-like and target-like fragments, this complex 
emission pattern can not be explained as a result of a 
simple statistical emission of the two excited participant 
nuclei. This fact suggests that it should be considered as 
a non-trivial phenomenon in which the interplay of the 
two-body interactions and the nuclear mean field plays a 
crucial role. 

There is not yet a uniform agreement about the possi- 
ble mechanisms responsible for this IMF overproduction 
at mid-rapidity. Recently, a statistical approach to the 
problem was developed [|ll| , and the claim that the inclu- 
sion of Coulombian interactions in the statistical picture 
could explain this behavior was made. On the contrary, 
some authors adopted a dynamical description to ana- 
lyze the same process. For example, a Landau- Vlasov 
semiclassical transport model was used to analyze the 
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presence Q and the temporal behavior of the par- 
ticle emission process that populates the mid-rapidity 
range. Neck instabilities, dynamical fluctuations, and 
particle production at mid-velocities were also studied 
using the quantum molecular dynamics (QMD) model 
in semi-peripheral ]l3t , and central heavy ion collisions 
[0. In addition, in a recent contribution, jl^], two mech- 
anisms of IMF production were suggested (neck, and 
fission- like surface emission), and two different models 
were introduced in order to describe them. 

The aim of this paper is to present a simple model that 
provides a description of the IMF production process at 
mid- velocities. In particular, we focus our attention on 
two mass asymmetric reactions of relatively small sys- 
tems, 58Ni-hC and ^^Ni+Au at 34.5 MeV/nucleon. The 
motivation to study these systems is two-fold. On one 
hand, it allows us to resolve the controversy about fast 
overlap emission processes and slow surface deformation 
breakup. On the other, our time-based analysis supports 
new experimental observations of alignment and proxim- 
ity effects in neck breakup of mass asymmetric collision 
partners fl^ . 

The reactions are studied within a quasi-classical ap- 
proach using molecular dynamics techniques. Out of var- 
ious methods used to study heavy ion reactions, the im- 
plemented MD model can describe changes of phase, hy- 
drodynamic flow, and non-equilibrium dynamics without 
adjustable parameters. All order nucleon-nucleon corre- 
lations to form fragments and neck-like structures are 
intrinsically incorporated in the dynamics. Of course, 
the price one has to pay using such a simplified quasi- 
classical description is that quantitative agreement with 
experiments can not be taken for granted. Nevertheless, 
despite its simplicity, we found it a great tool to probe 
new physical effects and to propose a meaningful scenario 
for different origins of MR particles (see [p^). 

Several kincmatical features of the experimental data 
are well reproduced. In addition, the availability of cor- 
relations of all orders at all times allowed us to study in 
detail the time evolution of the fragment formation pro- 
cess of IMF's and to identify different mechanisms of mid- 
rapidity particle production within a unified description. 
In this way, three time-scales could be differentiated: a 
violent pre-equilibrium stage, a delayed aligne emission 
stage, and a "statistical" evaporation stage. 

In section II, we briefiy describe the experimental 
setup. The description of the nuclear interactions model 
and the used fragment recognition algorithm are given 
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in sections III and IV respectively. The obtained results 
and their posterior analysis are included in sections V 
and VI. Finally, conclusions are drawn in section VII. 



II. EXPERIMENTAL DATA 

The two mass asymmetric reactions studied in this 
contribution were performed at the coupled Tandem 
and Super-Conducting Cyclotron accelerators of AECL 
at Chalk River. A beam of ^^Ni accelerated at 34.5 
MeV/nucleon bombarded alternatively a 2.4 mg/cm^ 
carbon target and a 2.7 mg/cm^ gold target. Charged 
particles issued from these reactions were detected in the 
CRL-Laval An array constituted by 144 detectors set in 
ten rings concentric to the beam axis and covering po- 
lar angles between 3.3° and 140°. See [|l^ and references 
therein for more information on detectors and energy cal- 
ibration. In the present work, we will consider events se- 
lected in the off-line analysis by a total detected charge of 
at least 24 and 34 (fully detected) units for the ^^Ni-|-Au, 
and ^^Ni-|-C systems respectively. 



III. THE MODEL 

In order to get information on the time formation and 
the kinematical characteristics of the intermediate veloc- 
ity material produced in these mass asymmetric heavy 
ion reactions, we relied on molecular dynamics simula- 
tions. The model used in this work was originally intro- 
duced by Lenk et al. to test the accuracy of the Vlasov- 
Nordheim approximation |18|. It is a classical molecular 
dynamics (CMD) model with a spherical two-body in- 
teraction potential in which Coulombian interactions arc 
taken into account, while the nuclear flavor is provided 
by the following interaction terms: 



Vnn=Vo[- 



Va\ 



(1) 
(2) 



Vnp and Vnn stand for neutron-proton and identical nu- 
cleons interactions respectively, and Vc — 5.4 fm is a 
cutoff radius. We used a set of parameters (M-set in 
Ref. [|l8|), that gives an equation of state with a com- 
pressibility around 250MeV, an equilibrium density of 
po = 0.16/TO-3, and E{po) = -IQMeVA. Due to the 
particular choice of parameters, this model provides a 
realistic nucleon-nucleon cross section and it has been 
used to reproduce reasonably well several features of ex- 
perimental data on heavy ion collisions (see ^ and 
references therein) 

The nuclei used in our computational experiments have 
been built with the right number of protons and neutrons 
while the corresponding ground-states were obtained by 



molecular dynamic techniques, i.e. by a cooling proce- 
dure that starts with a confined and rather excited nu- 
cleus, and ends with a self-contained state at reasonable 
binding energies. Before each collision a random relative 
projectile-target orientation was chosen, and afterwards 
the projectile's center of mass velocity was boosted to 
the desired energy. The set of equation of motion is inte- 
grated up to a final time oit — 2000 fm/c with a standard 
velo city- Ver let algorithm |^^, taking Unt = 0.02 fm/c as 
the integration time step, achieving an energy conserva- 
tion of 0.01%. Under this scheme, we have analyzed 8400 
cohisional events for the ^^Ni-|-^^C at 34.5 MeV/nucleon 
reaction, for a wide range of impact parameter values 
{h < 4.5 fm), and 3400 events for the ^SNi-h^^^Au at 34.5 
MeV/nucleon case (5fm< b < 8fm). 



IV. FRAGMENT RECOGNITION 

In order to explore the nucleon-nucleon correlations 
that give rise to the fragment formation process, we have 
adopted a definition of cluster compatible with the clus- 
terization algorithm known in the literature as the mini- 
mum spanning tree in energy space (MSTE) ||2^ . Under 
the MSTE scheme a given set of particles i,j,...,k be- 
longs to the same cluster Ca if: 



(3) 



where = V{i'ij) + (Pi ~ Pj)^/4Ai, and fi is the re- 
duced mass of the pair {i,j}. This clusterization method 
searches for correlated structures in q space (first term 
in the definition of e.y ) , with the additional possibility of 
avoiding certain particle incorporation to a given cluster 
regarding at the relative momenta of particle pairs (sec- 
ond term of e^). In this way the conformation of the 
MSTE partitions refiects certain degree of correlation in 
q-p space. Due to this feature, the MSTE becomes very 
useful in the recognition process of promptly emitted par- 
ticles during the pre-equilibrium stage of the reaction. 

The availability of correlations of all-orders at all times 
allows us to gain a lot of insight about the presence and 
origin of mid-velocity sources. For instance, right after 
the most violent stage of the collision, an MSTE clus- 
ter recognition step enables the tagging of particles as 
belonging to the quasi-projectile (QP), the quasi-target 
(QT), or as free particles (FP) not q-p correlated neither 
to the QT nor the QP at tagging time {ttag)- 

The value of ttag is determined on an event by event 
basis. Analyzing the temporal behavior of the size of the 
two biggest clusters, ttag is associated with the time at 
which the mass number of these fragments attains stabil- 
ity for the first time after the fragmentation process has 
started ||2^. Figure |l| illustrates this last point. It shows 
a typical ttag determination for a single ^^Ni -i-^^"^ An re- 
action with an impact parameter value of & = 6/to. 
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FIG. 1. Mass temporal evolution of the two biggest MSTE 
clusters for a typical ^^Ni Au reaction with an impact 
parameter value of 6 = 6/m. The dashed line signals the 
estimated ttag- 

It is then clear that ttag signals the entrance into a 
new and less violent deexcitation regime with particles 
QP(QT)-tagged being well correlated in configuration 
space right after the highly coUisional stage. 

In addition, the microscopic description allows to eas- 
ily calculate the time (te) at which every fragment has 
actually been emitted by simply tracing back the dy- 
namical evolution of the reaction. Since each asymp- 
totic detected cluster is emitted at different stages of the 
reaction (a cluster is considered as a whole when most 
of its nuclcon constituents get well correlated velocities 
and positions, differentiating themselves from other frag- 
ments) it is possible to extract useful information cor- 
relating this emission time characterization with other 
kinematical variables. 



V. RESULTS 

A. Reliability of the simulations 

In order to compare the numerical simulations with the 
experimental data the following prescription was adopted 
to estimate an observable (bexp) related to the reaction 
impact parameter (see Q, psj): 



MA 
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where nj^™ is the total parallel momentum of all charged 
particles in the forward velocity hemisphere of the center 
of mass (CM) reference frame, Ap, A-p and Pp™ are re- 
spectively the projectile mass number, target mass num- 
ber and projectile CM momentum. We took tq = 1.2/to. 
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FIG. 2. Normalized parallel velocity distribution, calcu- 
lated in the center of mass system for z > 3 particles, cor- 
responding to mid-peripheral and peripheral ^^Ni + C and 
^^Ni -\- Au reactions are shown in upper and lower pan- 
els respectively. Vertical dashed lines were drawn at nu- 
cleon-nucleon velocity value. Figures (a) and (c) show ex- 
perimental results, whereas (b) and (d) show simulated ones. 
In each panel the IMF contribution is shadowed, whereas QP 
tagged IMF's are stripped. See text for details. 

In panels (a) and (b) of figure ^ we show the normal- 
ized parallel velocity distribution (w[j'"), in the CM ref- 
erence frame, of the z > 3 collision charged products 
for '^^Ni-t-^^C peripheral and mid-peripheral [h^xp > 2.5 
fm) reactions. Experimental data is shown in the first 
panel and simulation results in the second one. The IMF 
contribution(3 < z < 7) was shadowed and quasi pro- 
jectile assigned IMF were stripped . (The Air detector 
coverage and the complete detected charge condition al- 
low us to establish, in this case, a meaningful comparison 
between the experimental data and unfiltered simulated 
reactions ). 

It can be seen that general features arc well repro- 
duced by the simulation. Aside from PLF and TLF ve- 
locity signatures, an intermediate bump peaked in be- 
tween Vnn and PLF velocity can be recognized for both, 
experimental and simulated data. Nevertheless, it seems 
that the CMD model overestimates slightly the amount 
of IMF's coming from the QP in the higher velocity re- 
gion of 0.03c. 

The corresponding normalized parallel velocity distri- 
butions for ^^Ni + Au mid-peripheral and peripheral re- 
actions (6 fm < bexp), are shown in panels (c) and (d) 
of figure |^. In this case, simulated data were filtered us- 
ing a software implementation of the detector. Despite 
the completely different entrance channel of this reaction, 
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and the different effect played by the detection thresh- 
olds of the apparatus, similar conclusions can be drawn 
for this reaction: general features of the distribution are 
well reproduced by the model, but an overproduction of 
IMF coming from the QP (centered at ~ 0.15c) can 
also be recognized in panel (d). 

We believe that this happens because our simplified 
model fails to accurately reproduce secondary deexcita- 
tion processes at a quantitative level. In particular, the 
relative amplitude of LCP's and IMF's production in the 
QP exit channel is not well reproduced because the quasi- 
classical approach can not adequately handle the relevant 
degrees of freedom involved in the slow nuclear deexcita- 
tion processes. 

Another point that is worth noting is that, as the dif- 
ference between the mid-rapidity bump (on the QP side) 
and Vnn velocities seems to be similar for both reactions 
(panels (a), and (c)), one can conclude that the deexci- 
tation process of the QP should be the same. Within our 
model, we can not rule out the possibility that in both 
reactions the QP proceed thru similar IMF production 
mechanisms. In particular, it seems possible that t he 
process of delayed asymmetric breakup (see Sec. VB 2 ) is 
there in both cases. However, along the present contri- 
bution, plausibility arguments will be given supporting 
the idea that this happens with different probabilities in 
each case. 

From the preceeding analysis, it can be concluded that, 
even if a quantitative description of the process seems to 
be questionable, the use of our CMD model can throw 
some light upon the problem of mid rapidity production. 
As will be shown, despite its simplicity, the model is flexi- 
ble enough to allow the presence of qualitatively different 
mechanisms of particle production that can be integrated 
in a unified description of the process. This can be used 
to settle a meaningful scenario within which experimen- 
tal data can be studied. 



B. CMD Description 

1. Source-origin analysis 

The microscopic description of the reaction allows us to 
associate two parameters to each asymptotic charged par- 
ticle: its emission time (ie), and its early physical origin, 
i.e if the particle belonged to the quasi-projectile (QP), 
quasi-target (QT) or was a prompt free particle (FP) at 



tag- 



In figure || we show the asymptotic distribution of t)|j™ 
calculated for every charged particle, for the two analyzed 
reactions. 
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FIG. 3. Normalized parallel velocity distribution, in the 
system center of mass, for every emitted charged particle in 
mid-peripheral ''**Ni-|-^^C (a), and ^**Ni-|-^^^Au (b) simulated 
reactions. The dashed vertical line corresponds to the nu- 
cleon-nucleon velocity value. The contribution of QP, FP, and 
QT-tagged particles are plotted striped, shaded and grided 
respectively. 

We considered, as in figure |^, mid-peripheral and pe- 
ripheral events. In the figures, the contribution of parti- 
cles that had been labelled as QT, FP, and QP at ttag, 
are distinctively shadowed. 

Figure || (a) corresponds to the ^^Ni-I-^'^C reaction 
(2.5fm< hf-xp < 6.0fm). It can be seen that the par- 
ticles tagged as FP were mainly emitted slightly above 
Vnn- Comparing with figure ^ (b), it can be noticed that 
this emission is almost entirely composed by LCP's, dy- 
namically expelled in the early stage of the reaction. The 
wide velocity distribution of the QT-tagged particles re- 
fiects the fact that the ^^C target is easily broken in the 
collision producing mainly z < 3 as byproducts. Finally 
it can also be seen, for the QP-tagged particles, an im- 
portant contribution of intermediate mass fragments and 
heavier particles, aside from the presence of LCP's (see 
figure H (b)). In this case, as was remarked earlier, not 
only a maximum centered at the projectile-like fragment 
(PLF) velocity can be seen, but a second peak can also 
be recognized towards the mid-velocity regions. 

Looking at figure (b) , a similar analysis can be per- 



formed for ^^Ni-h^^'^u collisions (9.0 fm< bexp < 12.0 
fm). In this case the FP parallel velocity distribution is 
centered at Vnn- Again this kind of emission is mainly 
composed of promptly emitted LCP's. The QP emission 
achieves its maximum at the PLF velocity, and presents 
a slight asymmetry towards the mid-velocity region. In- 
deed, this last feature is much more noticeable on the 
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side of the heavier partner of the reaction, where an ex- 
tra contribution towards the mid-rapidity, aside of the 
parallel velocity distribution centered around the target- 
like fragment (TLF) velocity, is easily recognizable. 



2. IMF's emission-time analysis 



As we mentioned at the end of section IV, in our CMD 
simulations, the emission time {tg) of asymptotic clusters 
can be easily estimated. At this point, it is important 
to remember that the tagged-origin classification reflects 
spatial and velocity correlations right after the compound 
nucleus starts to fragment, while signals the time when 
a given cluster attains mass stability. 

In figure ^ we used such characterization in order to 
plot 2-D velocity distributions, in the CM system, of 
IMF particles (3 < Z < 7) for ^^Ni-l-i^C mid-peripheral 
reactions (2fm < h^xp < 4fm). The figures show the 
velocity distribution in the perpendicular plane of the 



collision {vy 



plane), and a logarithmic scale was 



used for the contour levels. The contribution of QT, FP, 
and QP-tagged particles are shown in the first, second, 
and third column respectively, while particles emitted at 
times ie < 150 fm/c, 150 fm/c< t^. < 600 fm/c, and 
> 600 fm/c after the collision are shown in the first, 
second and third row respectively. In the figures, we in- 
cluded, as reference, the mean velocities (taken over all 
charged particles) of QP-tagged and QT-tagged particles 
(solid circles). 
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FIG. 4. Two-dimensional velocity distribution, for 
mid-peripheral ^**Ni-f ^^C simulated reactions. A cut perpen- 
dicular to the reaction plane is shown, i.e. {vy,Vz). The 
first, second, and third columns include data for QT, FP, 
and QP-tagged IMF's respectively. The first, second and 
third figure rows show the contribution of particles emitted 
at te < 150 fm/c, 150 fm/c< te < 600 fm/c, and te > 600 
fm/c after the beginning of the collision. 

It can be seen that the QT-tagged IMF particles are 
emitted very early in the evolution (Figure |4| (a)). They 
correspond mainly to small IMF's that come from the 
^^C breakup and, eventually, some dragged ^^Ni nucle- 
ons. On the other hand the FP-tagged particles (middle 
panels of figure that were already differentiated at te 
from the quasi-projectile and the quasi-target, seem to be 
emitted in a conical pattern for tf, < 150 fm/c reflecting 
the violence of the collision and the extremely asymmet- 
ric entrance channel. A second contribution of this kind 
of particles occurs a little bit later, 150 fm/c< te < 600 
fm/c, when excited promptly ejected FP-clusters emit a 
second wave of FP-fragments that populate the interme- 
diate velocity region, see Figure ^ (e). (It is important 
to keep in mind, though, that the main contribution of 
FP-tagged particles are not IMF's but LCP's) Finally, 
the QP-tagged particles, shown in the third column, ex- 
hibit a very interesting behavior. The maximum IMF 
emission of this kind of particles occurs at intermediate 
times, 150 fm/c< te < 600 fm/c, after the collision (Fig- 
ure ^ (f)). A Coulomb hole in the emission pattern, that 
reflects the interaction between the QP-residue and the 
emitted IMF clusters, can easily be recognized. This kind 
of pattern explains the second maximum observed in the 
parallel velocity distribution shown in Figure ||(a), and 
can be considered as the signature of a dynamically in- 
duced delayed QP breakup. Note also that after this kind 
of delayed emission, a subsequent isotropical pattern can 
be observed for longer times (Figure || (i)). 
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FIG. 5. Same as figure 4, but calculated for ^**Ni+^^''Au 
simulated reactions. 

The same analysis can be applied to the ^*Ni+^^^Au 
reaction (see Figure ||), where the same kind of behavior 
can be observed, for very short and very long emission 
times (first and third row respectively). Nevertheless, 
the emission pattern displayed by the system at interme- 
diate times, 150 fm/c< < 600 fm/c, after the colli- 
sion presents a special interest. One can observe that in 
this case, a delayed IMF emission pattern is recognizable 
mainly on the QT side of the reaction, while the QP re- 
lated IMF emission is nearly isotropic at this times, in 
contrast to figure |4|(f). 

We complete the characterization of the delayed IMF 
emission observed in both reactions with figure ^. This 
figure shows, for 150 fm/c< < 600 fm/c, the two- 
dimensional velocity distributions in the reaction plane 
(first column), and in the perpendicular plane (second 
column). The first and second rows of figures corre- 
spond to ^^Ni-I-^^C and ^*Ni-|-^'^''Au respectively. For the 
sake of clarity, we did not make any distinction between 
QP, QT or FP-tagged contributions in this case. It can 
be seen that the observed delayed IMF emission occurs 
mainly along the QP-QT direction in both reactions. A 
slight deviation of this alignment is observed due to the 
angular momentum gained by the emission sources after 
the collision. 
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FIG. 6. Two-dimensional velocity distributions for IMF 
particles emitted between te = 150fm/c and te = 600fm/c. 
Frames (a) and (b) correspond to mid-peripheral ^^Ni-I-'^'^C 
simulated reactions, while frames (c) and (d) to ®*Ni-|-^^^Au 
ones. The left column shows the reaction velocity plane, 
{vx,Vz), while the right column shows a perpendicular cut 

(Vy, Vz)- 



3. Role of Coulomb Interactions 

In this section we will analyze the specific role played 
by the Coulombian interactions in the detected delayed 
IMF emission. To that end we study the ^^Ni + Au 
reaction, comparing the behavior of the system when 
Coulomb interactions are not considered. 

In figure ^ parallel velocity distributions are shown 
when Coulomb is neglected (panels (a) and (b)) and when 
it is considered in the simulated evolutions (panels (c) 
and (d)). QT emitted particles and QP emitted ones are 
displayed on the first and second columns respectively. 
The contribution of z > 3 particles is stripped, whereas 
IMF contribution is cross-hatched. 
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VI. DISCUSSION 
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FIG. 7. Parallel velocity distribution of QT and QP emit- 
ted particles (panels (a),(c), and (b),(d)) in the simulated 
^^Ni+Au reaction. « > 3 particle contribution is stripped, 
whereas IMF contribution is cross-hatched. The upper panels 
show simulations without Colombian interaction. The lower 
ones include it. 

It can be seen from the figure that the QT emis- 
sion pattern is qualitatively affected by the presence of 
Coulombian interactions. In panel (a) the QT IMF pro- 
duction is negligeable. The few IMF present in that 
picture are barely attached to the quasi target , being 
realeased within a very short time after the projectile 
pass through. Without Coulomb instabilities, the dy- 
namically induced deformations are mostly reabsorbed 
by the TLF. 

This is not the case when electrostatic interactions are 
considered (panel (c)). Moreover, the overall QT emis- 
sion pattern in panel (a) is rather isotropic and have lost, 
up to a certain degree, the memory of the entrance chan- 
nel. On the contrary, the QT emission pattern when 
Coulombian interactions are included in the Hamiltonian 
(panel (c)), is severely affected by the collisional dynam- 
ics, giving place to what we called IMF aligned delayed 
emission. 

The QP emission pattern, on the other hand, shows 
the same qualitative behavior with or without Coulomb. 
That would imply that electrostatic interactions are not 
the key ingredient in the QP deexcitation process. 



All of these findings show the existence of two different 

scenarios for the origin of intermediate velocity particles. 
On one hand, there is an important prompt emission of 
LCP particles, at the highly collisional stage of the re- 
action. They come mainly from the overlap zone, as a 
consequence of nuclcon-nucleon collisions that promptly 
eject them from the 'bulk' nuclear mean field. 

On the other, a mechanism where a dynamical defor- 
mation of the heavier partner of the reaction develops can 
also be recognized. Eventually, this leads to the emission 
of the 'attached neck' followed by a Coulombian push 
(due to the proximity of the heavy source) that projects 
the emitted particles towards the intermediate velocity 
range. This second alternative occurs on a slower time 
scale and mainly involves the emission of IMF particles. 

This second type of emission can be considered as a dy- 
namically induced asymmetric fission. In fact, the study 
of two reactions with an inverse asymmetry in the en- 
trance channel allowed us to reinforce the idea that this 
kind of delayed aligned emission has its origin in a dy- 
namical induced shape instability that mostly affects the 
heavier partner of the reaction. More important deforma- 
tions can be sustained by the heaviest nucleus for enough 
time to allow asymmetric fission processes to develop. 

Coulomb effects were analyzed in the ^*Ni-|-Au reac- 
tion, and no qualitative changes were observed in the 
emission pattern of the QP, while major ones were re- 
ported for the QT. This fact reinforces the idea that 
Coulomb instabilities develops mostly on the side of the 
biggest partner, possibly not only because of the larger 
charge involved, but also because the dynamical defor- 
mation settled by the rescparation dynamics is more pro- 
nounced on the biggest nuclei side. 

When the two partners of the reaction reseparate, the 
biggest nucleus in an asymmetric reaction uses more of 
its surface nucleons than do the smallest one to estab- 
lish the necklike bond between the two poles. Therefore, 
the induced shape deformation on the biggest nucleus 
after rescparation should be greater (the smallest nu- 
cleus has relatively not enough surface to reabsorb the 
neck). The neck either proceed through multiple neck 
ruptures before reabsortion, breaking up in the reported 
aligned asymmetric way, or is effectively reabsorbed by 
the biggest nucleus. Moreover, if we assume an equal en- 
ergy sharing picture, valid for peripheral events, the Ni 
nucleus would be much hotter than the Au one and could 
breakup by other processes before shape deformations 
develops. Even if some fragments could still be emitted 
by dynamical deformations in this case, the associated 
partial cross sections would be very hard to isolate. It 
is worth noting, though, that even if the presented pic- 
ture favors the occurence of IMF delayed emission on the 
biggest nucleus side, nothing forbids the phenomenon to 
appear on a less important scale on the smallest nucleus 
side. 
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VII. CONCLUSIONS 



In this paper, we have studied two asymmetric heavy 
ion reactions using the description scheme of a classical 
molecular dynamics model. We found that calculations 
are consistent with experiment, and that processes de- 
scribed by the model could be relevant to the experimen- 
tal case. 

We have made extensive use of the availability of the 
microscopic correlations at all times provided by this ap- 
proach. Adopting the MSTE definition of fragment, the 
detailed knowledge of the dynamics allowed us to classify 
the complete set of asymptotic charged particles accord- 
ing to its physical situation in phase space, right after 
the collision (QP, QT or FP classification). Correlating 
this information with the emission timescale associated 
to each fragment, a clear distinction between different 
modes of particle emission in the mid-rapidity range was 
established. A fast emission process that mainly involve 
LCP, and a delayed aligned emission of IMF's, mainly 
coming from the heaviest partner of the reaction, were 
identified. 

The delayed aligned emission pattern reported in this 
paper could possibly be compared to the statistical emis- 
sion of reference ||ll|, where a nucleus deexcite in the 
presence of the Coulomb field of a secondary source. Re- 
sults of Ref. were obtained by simulating a Au-|-Au 
reaction at 35 MeV/nucleon, therefore in the presence 
of strong Coulomb fields. According to our calculations, 
the same phenomenon of delayed aligned emission can 
be traced even in the ^*iVi-|-C reaction, despite the small 
size of the TLF. Experimental evidences for this effect 
were also reported in jl^]. This could indicate that dy- 
namical effects can not be disregarded in the description 
of the process for light systems. 

All of these findings set a reasonable framework in 
which recent experimental observations of particle align- 
ment and proximity effects of the heavier partner can 
be understood Moreover, they should be taken into 
account in the interpretation of other experimental obser- 
vations. For example, in what concerns the geometrical 
and shape constraints that the observed delayed emission 
of IMF's impose to chemical or thermal equilibration pro- 
cesses in nuclear systems. 
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